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Abstract. We consider the model Poisson problem —Au = / £ Q, u = g on 
dfl, where H is a bounded polyhedral domain in R n . The objective of the paper 
is twofold. The first objective is to review the well posedness and the regularity 
of our model problem using appropriate weighted spaces for the data and the 
solution. We use these results to derive the domain of the Laplace operator 
with zero boundary conditions on a concave domain, which seems not to have 
been fully investigated before. We also mention some extensions of our results 
to interface problems for the Elasticity equation. The second objective is to 
illustrate how anisotropic weighted regularity results for the Laplace operator 
in 3D are used in designing efficient finite element discretizations of elliptic 
boundary value problems, with the focus on the efficient discretization of the 
Poisson problem on polyhedral domains in R 3 , following Numer. Fund. Anal. 
Optim., 28(7-8):775-824, 2007. The anisotropic weighted regularity results 
described and used in the second part of the paper are a consequence of the 
well-posedness results in (isotropically) weighted Sobolev spaces described in 
the first part of the paper. The paper is based on the talk by the last named 
author at the Congress of Romanian Mathematicians, Brasov 2011, and is 
largely a survey paper. 



Introduction 

Let CI C R™ be an open, bounded set. Consider the boundary value problem 

f Au = f in Cl 

[u\dn=g, on £2, 

denned on a bounded domain Cl c R™, where A is the Laplacian A = Ylt=i^l- 
When dCl is smooth, it is well known that this Poisson problem has a unique solution 
u G H m+1 (Cl) for any / G J?™- 1 (ft) and g G H m+1 / 2 (dCl) [33 H3 [53J . Moreover, 
u depends continuously on / and g. This result is the classical well-posedness of 
the Poisson problem on smooth domains. 

On the other hand, when Cl is not smooth, it is also well known [331 IM1 E3 139 > 40 
that there exists s = sq such that u G H s (Cl) for any s < sq, but u g" H sn (Cl) in 
general, even if / and g are smooth functions defined in a neighborhood of CI. For 
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instance, if Q is a polygonal domain in two dimensions, then sq = 1 + 7r / olmax i 
where Umax is the largest interior angle of f2 [59]. See also Wahlbin's paper [53] . 

In view of applications to the Finite Element Method, we restrict out attention 
to domains with polyhedral structure. These are natural non-convex generalizations 
of classical n-dimensional polyhedra that allow for curved boundaries, cracks (i.e., 
internal faces), and non-smooth boundary points touching a smooth part of the 
boundary. We refer to [T3] for a precise formulation. 

There exist a very large number of papers devoted to boundary value problems 
on non-smooth domains. While it is impossible to mention them all, let us at least 
mention the papers of Arnold, Scott, and Vogelius [TU], Babuska and Guo [55] . 
Bacutja, Bramble, and Xu [19], Jerison and Kenig [37], Kondratiev [39], Kozlov, 
Mazya, and Rossmann [30], Mitrea and Taylor [37], Rossmann [5T], Verchota [53], 
and many others. Other results specific to numerical methods for polyhedral do- 
mains are contained in the papers of Apel and Dobrowolski [5J, Costabel, Dauge, 
and Nicaise [25], Costabel, Dauge, and Schwab [23] ) Dauge [23], Demkowicz, Monk, 
Schwab, and Vardapetyan [25], Elschner [26], Guo and Schwab [35], and many oth- 
ers. Further results and references can be found in the aforementioned papers, as 
well as in the the monographs of Grisvard [52] as well as the recent book [39"] . 
Regularity for polyhedral domains is useful in designing fast solvers for numerical 
methods [7J[T7]. See also [51 El [12] [TSJ [THl Ell ESI 11H1 Sil SI] for more applications 
of these techniques to other types of Partial Differential Equations and numerical 
methods. 

In this paper, we shall review the results from [T31 [T3J [HE 133] , and [3S], which 
make use of the natural stratified space structure on f2. This leads, by successive 
conformal changes of the metric, to a metric for which the smooth part of f2 becomes 
a smooth manifold with boundary whose double is complete. The resulting Sobolev 
spaces defined by the new metric will lead to spaces on which the Poisson problem 
is well-posed. 

We restrict for simplicity to consider only the Laplace operator in ([I]). However, 
all theoretical results presented here extend to scalar, strongly elliptic, linear opera- 
tors P with sufficiently regular coefficients, and even to elliptic systems, such as the 
system of anisotropic elasticity [46] - Furthermore, we can also treat transmission 
problems, for which the coefficients of P are allowed to jump across piecewise- 
smooth hypersurfaces, representing interfaces, under some additional conditions 



[131 132] • We briefly discuss these extensions in Subsection 1.3 



For the discretization on polyhedral domains, we build discrete spaces Sk C 
Hq(Q) and Galerkin finite element projections Uk S Sk that approximate the so- 
lution u of Equation for / e iJ m_1 (f2) arbitrary. We prove that, by using 
certain spaces of continuous, piecewise polynomials of degree to, the sequence Sk 
achieves quasi-optimal rates of convergence. More precisely we prove the existence 
of a constant C > 0, independent of k and /, such that 

(2) ||u - u k \\ m(n) < Cdim(5 fc )- m/n ||/||ir»-i(n), «* G S k , 

where n = 2 or n = 3 is the dimension of our polyhedral domain. 

The contents of the paper are as follows. In the first section we review well- 
posedness results in weighted Sobolev spaces on polyhedra domains. These weighted 
spaces are sometimes called the Babuska-Kondratiev spaces. These results are not 
sufficient for our applications to the Finite Element Method in three dimensions, 
so in the second section we review some additional anisotropic regularity results. 
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These results are used in the third section to construct a sequence of meshes that 
yields /i m -quasi-optimal rates of convergence in three dimensions. Finally, in the 
last section we discuss some of the main ingredients that enter in the proof and 
which are of independent interest. These include the Hardy-Poincare inequality 
(which guarantees the coercivity of our problem) and a description of the weighted 
Sobolev (or Babuska-Kondratiev) spaces JC™(Sl), which are the natural spaces for 
our well-posedness results, as the usual Sobolev spaces for a modified metric on f2, 
which nevertheless is conformally equivalent to the old one. 

Acknowledgements: The authors would like to thank the organizers of the Inter- 
national Congress of Romanian mathematicians, Brasov 2011, where these results 
were presented. A. M. would also like to acknowledge the support of the Mathe- 
matics Department at Cornell University, where she is currently a Michler Fellow. 

1. Well posedness in isotropic weighted Sobolev spaces 

Using the standard notation for partial derivatives, namely dj = J^r and d a = 
d" 1 . . . d" n , for any multi-index a = {ct\ , . . . , a n ) S , we denote the usual Sobolev 
spaces on an open set V C R" by 

H m (V) = {u : V -> C, d a u E L 2 (V), \a\ < m}. 

As mentioned in the introduction, the solution of our model Poisson problem ([!]) 
has only limited regularity in the spaces H m (Q). The situation changes for the 
better if one considers weighted Sobolev spaces, though. To define the weighted 
analogues of these spaces, we need to introduce the notion of singular boundary 
points of the domain C W 1 . 

1.1. Weighted Sobolev spaces. Let d smg ft C d£l be the set of singular (or non- 
smooth) boundary points of fi, that is, the set of points p E dVl such <9f2 is not 
smooth in a neighborhood of p. In case we consider mixed boundary conditions, the 
set of singular points includes also the set of points where the boudary conditions 
change. If, furthermore, interfaces are considered, the set of singular points contains 
the set of singular points of the interface, as well as the set of points where the 
interface touches the boundary. We will denote by rn (x) the distance from a point 
iGlJto the set 9 s i ng f2 and agree to take rn = 1 if there are no such points, i.e., if 
dil is smooth. For /i E Z + and a E R, we define the weighted Sobolev spaces 

(3) JC£(fi) = {uE Ll c (Q), r^- a d a u E L 2 (tt), for all \a\ < M }, 

which we endow with the induced Hilbert space norm. We note that for n = 3 
for example and a polyhedral domain in ]R 3 , we have that rn(x) is the distance 
to the skeleton comprising the union of the closed edges of dil. Recently, general 
spaces of this kind were studied by H. Amann [HQ]. 

Similar weighted Sobolev spaces are associated to the faces of ft. By a face, we 
mean the connected components of the boundary dfl after the set of singular points 
is removed. For example for n — 3, we define 

= {(u F ), r^~ a d a u F E L 2 (F) }, 

where \a\ < m and F ranges through the set of faces of dfl. For s E R+, we define 
the space K-KdCl) by standard interpolation. 
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1.2. Well-posedness for the Poisson problem on n-dimensional polyhedral 
domains. The following result is proved in (T3]. For simplicity, we shall assume 
that Q has no cracks and that there are no vertices that touch the boundary. (That 
is, we shall consider only domains fl that coincide with the interior of their closure 
Q.) 

Theorem 1.1. Let £1 C K™, be a bounded, curvilinear polyhedral domain and m £ 
Z_|_. Then there exists t]q > such that A(w) = (Am, u|sq) defines an isomorphism 

/or |a| < 770. If m = 0, </ie solution u corresponding to the data (/, 0) G 
/C^-^fi) /C*^ 2 ^ (<9f2) is aLso £/ie solution of the associated variational problem. 

This theorem amounts to the well-posedness of the boundary value ([TJ on n- 
dimensional polyhedral domains. For n — 2 (in which case is a polygonal domain), 
this result is due to Kondraticv 391, in which case 770 = — - — , where olm ax is the 

I S> O.MAX ivitxj\ 

measure in radians of the maximum angle of f2. For n = 3, this result was proved 
in []3]. For later applications, we shall need the following result. 

Theorem 1.2. The results of Theorem \l . 1\ remain true for infinite angles in two di- 
mensions, infinite polyhedral cones in three dimensions, and infinite dihedral angles 
in three dimensions. 

The proof of this theorem proceeds along the lines of the proof of Theorem |1.1| in 
[2] or [X3] . A first difference to remark is that the Hardy-Poincare inequality does 
not hold for the whole domain. Then the "desingularization" S(57) has to involve 
the directions at infinity also in the case of an angle or a cone. In case of a dihedral 
angle D a — {0 < 9 < a}, in cylindrical coordinates (r, 6, z), one has to consider 
the also the two point compactification of the edge. In particular, 

(4) £(D a ) = [Q, a] x [0,oo] x [-00,00]. 

1.3. Extensions. Theorem |1.1| above was extended in several ways. First, the 
proof applies with almost no change if mixed boundary value problems are con- 
sidered, provided that no adjacent faces are endowed with Neumann boundary 
conditions. We do allow, however, different boundary conditions on the same face. 
We treat the points where the boundary conditions change similarly to the non- 
smooth boundary points, as solutions exhibit a similar singular behavior in this 
case. 

As already mentioned in the Introduction, we can more generally consider a 
general uniformly strongly elliptic differential operator of the form 

(5) Lu = — di(aijdju) + cu, with c > 0. 

ij 

(Recall that L is uniformly strongly elliptic if, and only if, there exists C > such 
that dijtitj >C^tl for all (U) £ R n .) 

We can also include certain transmission or interface problems. More precisely, 
we now assume that our domain fl can be written as a union of curvilinear poly- 
hedral domains Vlj with disjoint interiors: £1 = uj^fl,. Let T := U^ =1 dftj \ Oil 
be the interface. We assume that T is smooth and assume further that no adjacent 
faces of the O^-'s are endowed with Neumann boundary conditions. We do allow 
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r to touch the boundary of f2. We can then extend the result of Theorem |1.1| by 
using instead the broken weighted Sobolev spaces /C™(f2), defined by 

(6) := ffif =1 /C™(^). 

We observe that, if there is no interface, = K™(p). We let 8d& be the 

part of the boundary with Dirichlet boundary conditions, which we assume to 
be a closed subset of the boundary, and let c^fi ■— dft \ Soil. We denote the 
outer normal vector to 57, which is defined a.e. on dfl, by v, and the conormal 
derivative associated to the operator L by — Via 13 ' {x)dj ;. Let L{u) — 

(Lu,u\d D Q, D^u\g N n). Our most general result in n dimensions states that for 
m > 1 L is an isomorphism (see |13]V 

(7) L : v a -j. it™-? (o) ffi /c™ 4 ^ 2 (fen) e /c™" 1 // (a^n) , 

where 

(8) V a := { u € n ^ +1 (fi), u+ = u", D^+u = Df-u on T }, 

and the subscript ± refers to non-tangential limits to each side of the interface. 
The conormal derivative is defined in the sense of the trace a.e. on dSl. 

Let us mention that the interface L will separate different faces where it touches 
the boundary, and hence we assume that these faces are not both endowed with 
Neumann boundary conditions. 

For elasticity with mixed boundary conditions, a similar result is obtained by 
Mazzucato and Nistor in 0B]. The results in [JS] also extend to interface problems 
under the same assumptions (no adjacent faces with Neumann boundary conditions 
and a smooth interface) using the methods as in [13] and in |46j . More precisely, 
we use Korn's inequality to obtain local regularity results (no weighted spaces). 
This applies, in particular, to interface problems. There the additional regularity 
is proved as for the additional regularity at the boundary for smooth domains. 
See [50] for a proof of the additional regularity at the boundary for systems that 
extends to interface problems. Once one has the local regularity results, the global 
regularity results in weighted spaces is proved as in [46] using suitable partitions 
of unity. The solvability in H 1 is an immediate consequence of Korn's inequality 
and of the Hardy-Poincare inequality. Combining regularity with solvability in H 1 
yields solvability in higher weighted Sobolev spaces /C^j^ 1 (SI) . 

Other regularity results go toward analytic regularity using countably normed 
spaces as in the work of Babuska-Guo [3HIS3], and Costabel, Dauge, and Nicaise 
|22) . See the Introduction for more references. It would be interesting to extend 
these results to the de Rham complex [U |9] . 

1.3.1. Adjacent Neumann faces and non-smooth interfaces in 2D. The assumption 
that no vertex P be the common point of two adjacent faces with Neumann bound- 
ary conditions or the assumption that T be smooth at any interface point P are 
both equivalent to the fact that the function constant equal to one not be a singular 
function at that singular point P. This assumption is necessary, because, if it is 
not satisfied, the relevant operator L is not even Fredholm for the value a = and 
it is also not invertible for any a G K. However, this assumption is not realistic 
in practice and, it turns out, not even necessary for designing graded meshes that 
yield quasi-optimal rates of convergence [53] . 
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To obtain a well-posedness result for interface problems in 2D, we can proceed 
as follows [13] ■ Let xp be a smooth function that is equal to 1 near each singular 
point P that is either a point where we have Neumann-Neumann conditions or a 
non-smooth interface point satisfying respectively Neumann or periodic boundary 
conditions on the sides at P. This includes points P that belong to more than two 
of the subdomains VLj (so called multiple junction points). We assume the x's have 
disjoint supports. Let W s be the linear span of the functions xp- The choice of 
boundary conditions or the introduction of additional singular points to a polygonal 
domain define a polygonal structure on 17, see |43j for details. 

Theorem 1.3. Let il be a domain with a polygonal structure. The there exists 
n > such that, for all any < a < rj and m G 1 + , the map 

with T> a given in is an isomorphism. 

The proof requires the calculation of the index of the operator L acting on 
^ ^-a+i(^)- Note that our result is not valid for a = 0. We expect a 
similar result in 3D. 



Theorem 1.3 can be used to justify the construction of a sequence of meshes (in 
2D) that yields quasi-optimal h m rates of convergence for transmission problems 
with non-smooth interfaces (and even with multiple junctions) and problems with 
adjacent Neumann-Neumann corners in 2D. See [JS] for additional issues related 
to the regularity and numerical methods for interface problems. We notice that 
the resulting sequence of meshes is the same for all 2D problems on polygonal 
domains (with or without interfaces or Neumann-Neumann corners), although the 
theoretical PDE result (or a priori estimates) are different in these two cases. 

1.4. The domain of A on concave polygons. Let us mention that the method 



used to obtain Theorem 1.3 can be used to describe the domain 2?(A) of the 
Fricdrichs extension of the Laplace operator on 11 with zero boundary conditions. 
First of all, the form associated to A, namely B(u,v) = (Vu, Vw), u,v zero on 
the boundary, defines the so called energy norm: \u\jji(Q) = _B(u, it) 1 / 2 . The com- 
pletion of C£°(f2) in the energy norm is i/o(ll). The proofs in [301 El show that 
Hq^CI) — K.{(Q)r\{u\gn — 0}, with equivalent norms. The domain of the Friedrichs 
extension of the Laplacian A is then 

(9) P(A) = {u G #d(fi), Au G L 2 {Q)}. 

If 17 is convex, then it is known that 2?(A) = if 2 (17) n Hq(VL). This is however not 
true if 17 is concave. To describe 2?(17) in the case when f2 is concave, let us notice 
that the map 

(10) A : /c|(si) n H^(n) -> L 2 (Vt) 

is Fredholm and its index is the number of re-entrant corners by |39) . Let P be 
such a re-entrant corner with angle ctp > n. Also, let (r, 9) be polar coordinates at 
P and consider the function (f>p — r 7T / ap sin(7r(9 / ctp)xp , where xp is the function 
considered in Theorem |1.3| Let V s be the space of linear combinations of the 
functions <pp, with P a re-entrant corner. Then one has that 

(11) A-.ic 2 2 (n)nH^(n) + v s ->L 2 (n) 

has index zero, is injective, and hence bijective. This proves the following result. 
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Theorem 1.4. The domain of the Friedrichs extension of the Laplace operator with 
zero boundary conditions on a polygon SI C K 2 is 

V(A) = /cl(si) n H^(Cl) + v s . 

A similar description is available for other types of boundary conditions. This 
result immediately leads to a maximal regularity result for the heat equation on 
polygonal domains. 

See also |3T] and [3D] for related results on Friedrichs extensions of second order 
elliptic operators on manifolds with conical points. 

2. Anisotropic weighted Sobolev spaces and regularity 

The well-posedness result of the previous section are not enough to establish 
quasi-optimal rates of convergence in 3D. We need additional regularity along the 
edges, as follows. Let u be the solution of problem with / e i? m_1 (n) and 
g = 0. We observe that this assumption is stronger than assuming that / is in 
a weighted Sobolev space of the form K,™~l for \a\ small. We will need to take 
advantage of this additional regularity of /, which leads to improved regularity 
for u along the edges. We encode this additional regularity by introducing new 
anisotropically weighted spaces. 

We assume first that the domain SI is a dihedral angle with axis along the in- 
coordinate axis, D a = {0 < 9 < a}, using cylindrical coordinates (r, 9, z). We 
further assume that / G H m ^ 1 (D a ). Then / G /C^ 1 (D a ), and hence 

(12) u€KZtf(D a ) 



for positive and small enough a, by Theorem 1.2 Hence, d z u £ K™(D a ). However, 
we also have Ad z u = d z Au = d z f G H m ~ 2 (£l). Then, using Theorem 1 1 . 1 1 which 
extends to this setting, we also obtain that 

(is) a,«6^i(fi), 



a better estimate than in Equation (12 1. These calculations suggest that we intro- 
duce a scale of spaces P™, m G Z+, as follows: 

T>\{D a ) :=K\{D a ), 

V™(D a ) := {u G K?{D a ), 8 z u G V™-\D a )}. 

The spaces T>\ are thus independent of a. 

We assume next that the domain SI is a cone C centered at the origin. Wc let 
p(x) = \x\, the distance from x to the origin, and define 

V\{C) := p^KliC) = {p^v, v G K.\{C)}. 

To introduce the spaces I?™(C) for m > 2, we shall need to consider the vector field 
pd p := xd x + yd y + zd z , which is the infinitesimal generator of dilations centered 
at the vertex of the cone. We then define by induction 

V™{C) := {u e JCT(C), pd p (u) G V™-\C)}, m > 2. 

For a general bounded polyhedral domain SI, we define the anisotropic weighted 
Sobolev spaces 2?™(S1) by localizing around vertices and edges, using as models 
cones and dihedral angles respectively, such that away from the edges these spaces 
coincide with the usual Sobolev spaces H m . Then, we have the following regularity 
result [T5] : 
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Theorem 2.1. Let f G H m (Ct), with in > 1. Then there exists ryo,a > such 
that the Poisson problem ([I]) with g — has a unique solution u G (f2) for 

any < a < r\ = r/fi and 

IMIz^+Hn) < ^,a||/|U"»-i(n)- 
See [H 122 HI EE] for related results. 



3. Quasi-optimal h m -MESn refinement 

We describe in this section a strategy to obtain quasi-optimal ft, m -mesh refine- 
ment. We follow [15] . from where the pictures are taken. The theoretical jus- 
tification of this construction is based on the anisotropic regularity result of the 
previous section, Theorem |2.1| Given a bounded polyhedral domain and a pa- 
rameter k G (0,1/2], we will provide a sequence T n of decompositions of Q into 
finitely many tetrahedra, such that, if S n is the finite element space of continuous, 
piecewise polynomials on T n , then is the Lagrange interpolant of u of order m, ur, n , 
has "quasi-optimal" approximability properties. The result can be formulated as 
follows: 

Theorem 3.1. Let a G (0, 1/2] and < k < 2~ m l a . Then there exists a sequence 
of meshes T n and a constant C > such that, for the corresponding sequence of 
finite element spaces S n , we have 



ui^Ih 1 ^) < C2 



-km I 



for any u € T>™_^~{£1), u\gn = 0, and for any k G Z + . 

Theorem ([2| is now a direct consequence of Theorem 3.1 



3.1. Refinement Strategy. Our refinement strategy will first generate a sequence 
of decompositions 7^ of f2 in tetrahedra and triangular prisms, while our meshes 
T n will be obtained by further dividing each prism in T„ into three tetrahedra. We 
now explain how the divisions 7^ are defined inductively, T^+i being a refinement 
of T n in which each edge is divided into two (possibly unequal) edges. 

To define the way each edge of 7^ is divided, we need to establish a hierarchy 
of the nodes of 7^. Therefore, given a point P £ fl, we shall say that P is of type 
V if it is a vertex of ft; we shall say that P is of type E if it is on an open edge 
of ft. Otherwise, we shall say that it is of type S (that is, a "smooth" point). The 
type of a point depends only on and not on any partition or meshing. The initial 
tetrahedralization will consist of edges of type VE, VS, ES, EE:=E 2 , and S 2 . 
We shall assume that our initial decomposition and initial tetrahedralization were 
defined so that no edges of type V 2 :— VV are present. The points of type V 
will be regarded as more singular than the points of type E, and the points of type 
E will be regarded as more singular than the points of type S. All the resulting 
triangles will hence be of one of the types VES, VSS, ESS. Let us notice that once 
our initial refinement is fine enough, the edges of our domain will be decomposed 
into segments of type VE and EE, and the segments of type EE will be containted 
in triangular prisms. Therefore, we can assume that there are no triangles of type 
EES. 

Our refinement procedure depends on the choice of a constant k g (0, 2 _m / a ), 



where a > is as in Theorem 2.1 and n < 1/2. We can improve our construction 



by considering different values of n associated to different vertices or edges. This 
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generalization can easily be carry out by the reader. See |43j for example. Let 
AB be a generic edge in the decompositions T n . Then, as part of the T n +i, this 
edge will be decomposed in two segments, AC and CA, such that \AC\ — k\AB\ if 
A is more singular than B (i.e., if AB is of type VE, VS, or ES). Except when 
K = 1/2, C will be closer to the more singular point. This procedure is as in [20] . 
See Figure |XT| 



A more singular than B 
\AC\ = k\AB\, k = 1/4 



A and B equally singular 
\AC\ = \AB\ 



Figure 3.1. Edge decomposition 



The above strategy to refine the edges induces a natural strategy for refining 
the triangular faces. If ABC is a triangle in the decomposition 7^, then in 7^ +1 , 
the triangle ABC will be divided into four other triangles, according with the edge 
strategy. The decomposition of triangles of type S 3 is obtained for k = 1/2. The 
type VSS triangle decomposition is described in Figure 3.2 (a). In the case when 
ABC is of type VES, however, we shall use a different construction. Namely, in 
this case we remove the newly introduced segment that is opposite to B, see Figure 
3.2 (b), and divide ABC into two triangles and a quadrilateral. The resulting 
quadrilateral will belong to a prism in 7^' +1 . 




A' c 

(a) A of type V or E 
(b) B and C of type S, \A'B\ = \A'C\ 
\AC'\ = k\AB\, \AB'\ = K \AC\ 



VER decomposition: ZE = 90° 
\VC'\ = k\VE\, \VB'\ = k\VR\ 
\EA'\ = k\ER\, A'C was removed 



FIGURE 3.2. Triangle decomposition, k = 1/4 



3.2. Divisions in tetrahedra and prisms. We now describe the construction of 
the sequence of the decompositions 7^ for n > 0. The required sequence of meshes 
T n will be defined by dividing all the prisms in 7^ into tetrahedra. For the first 
level of semi-uniform refinement of a prism, more details are presented in [15] . 

We start with an initial division 7~q of £1 in straight triangular prisms and tetra- 
hedra of types VESS and VS 3 , having a vertex in common with Q, and an interior 
region Ao. See Figure 3.2 (a), where we have assumed that our domain Q is a 
tetrahedron. For each of the prisms we choose a diagonal (called mark) which will 
be used to uniquely define a partition of the triangular prism into tetrahedra. We 
then divide the interior region A into tetrahedra that will match the marks. Also, 
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(a) Initial decomposition. (b) Marking a prism: BC — mark, 

AA' || BB' || CC _L ABC and A'B'C 



Figure 3.3. The initial decomposition 7q of fl. 



we assume that the marks on adjacent prisms are compatible, so that the resulting 
meshes are conforming. We can further assume that some of the edge points (as in 



Figure 3.2 (b)) have been moved along the edges so that the prisms become straight 
triangular prisms i.e., the edges are perpendicular to the bases. 

The decompositions 7^ are then obtained by induction following the Steps 1 
through 3 explained next. We assume that the decomposition 7^ was defined and 
we proceed to define the decomposition T T [ +1 . 

Step 1. The tetrahedra of type S 4 are refined uniformly by dividing along the 
planes given by Xi + Xj = k/2 n , I < k < 2", where Xj are affine barycentric coordi- 
nates. This refinement is compatible with the already defined refinement procedure 
for the faces. See Figure [3~2| (a) for n = 1. 

Step 2. We perform semi-uniform refinement for prisms in our decomposition 7^ 
(all these prisms will have an edge in common with f2). This procedure is shown in 



Figure 3.2 (b) 





(a) First level of uniform refinement (b) First level of semi-uniform refinement 

of a prism, CD = mark 



FIGURE 3.4. First refinement T{. 



Step 3. We perform non-uniform refinement for the tetrahedra of type VS 3 and 
VESS. More precisely, we divide a tetrahedron of type VS 3 into 12 tetrahedra as 
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in the uniform strategy, with the edges through the vertex of type V divided in the 
ratio k. We thus obtain one tetrahedron of type VS 3 and 11 tetrahedra of type 
S 4 . (At the next step, which yields 7^ +2 we iterate this procedure for the small 
tetrahedron of type VS 3 , while the tetrahedra of type S are divided uniformly.) 
See Figure 3.2 (a). On the other hand, a tetrahedron of type VESS will be divided 
it into 6 tetrahedra of type S 4 , one tetrahedron of type VS 3 , and a triangular 



prism. The vertex of type E of will belong only to the prism. See Figure [372] (b). 
This refinement is compatible with the earlier refinement of the faces. 




The description of our refinement procedure is now complete. 



3.3. Intrinsic local refinement. We see that one of the main features of our 
refinement is that each edge, each triangle, and each quadrilateral that appears in 
a tetrahedron or prism in the decomposition 7^ is divided in the decomposition 
T.,' l+ i in an intrinsic way that depends only on the type of the vertices of that edge, 
triangle, or quadrilateral. In particular, the way that a face in 7^ is divided to yield 
Tn+i does not depend on the type of the other vertices of the tetrahedron or prism 
to which it belongs. This ensures that tetrahedralization T n +i, which is obtained 
from 7^ +1 by dividing each prism in three tetrahedra, is a conforming mesh. 



4. Hardy-Poincare inequality and regularity: a glimpse at the 

PROOFS 

There are two main ingredients for the proofs of the well-posedness results stated 
in the first section. One is the Hardy-Poincare inequality, which yields solvability 
(more precisely well-posedness) in the H^-type spaces and the second one is a 
regularity result, which allows us then to obtain well-posedness in higher regularity 
spaces. A third, more technical ingredient, is to describe the trace spaces at the 
boundary. For this, we use the same ideas as the ones used in the proof of regularity. 
We now discuss these ingredients. 
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4.1. The Hardy-Poincare inequality. Let us denote by tq(x) the distance from 
x to the set of singular points in the boundary of ft. Recall that these singular 
points consist not just of the edge points, but also of the points where the bound- 
ary conditions change and the points where the interface touches the boundary. 
The following inequality is then proved by induction [13] (see [14] for the three 
dimensional case, the two dimensional case was well known, see |49j for example). 

Proposition 4.1. Let ft be a polyhedral domain in W 1 . We assume that either 
ft is bounded, or that it is a cone or a dihedral angle. Let us assume that the 
Neumann part of the boundary d^ft := dfl \ doft contains no adjacent faces of ft. 
Then there exists a constant Cn > 0, which depends only on ft and the choice of 
boundary conditions such that the following Hardy-Poincaree inequality holds: 

[ ^dx <C n [ \\/u\ 2 dx 

for any function u G i? x (f2) that is zero on do ft- 

Let us assume that ft is bounded. A simple consequence of the Hardy-Poincare 
inequality of Proposition |4.1| is that the spaces 

(14) Hj } (ft) = {ueH 1 (fl),u = Oond D fl} 
and 

(15) JC\{fl) n{u£ Hl oc (Cl),u = on d D ft} 

are the same and their respective norms are equivalent. Neither this result nor the 
Hardy-Poincare inequality are true if there exist two adjacent faces with Neumann 
boundary conditions. This is the reason we needed a different approach in Section]!] 

4.2. Sobolev spaces and regularity. Our definition of weighted Sobolev spaces, 
Equation ([3]), is elementary. However, for the purpose of establishing the needed 
properties of these spaces, it is convenient to identify them with the usual Sobolev 
spaces associated to a different metric on ft. 

To this end, let us recall from |13j that a stratified curvilinear polyhedral domain 
ft is an open subset of a Riemannian manifold (M , g) of dimension d together with 
a stratification of 

(16) ft = n {d) D fl^ 1 ^ D . . . D D fi<°>. 

We then define stratified curvilinear polyhedral domains by induction as follows. 

For d = 0, ft is just a finite set of points. For d = 1, ft is a finite set of intervals. 

The stratum So for d — 1 will contain all the boundary points of the intervals, but 

may contain also other points. For d > 1, we require our domain Q to satisfy the 

following conditions: for every point p € dfl, there exist a neighborhood V p C M 

such that, if p G fl^ \ f^ -1 * 1 , I = 1, .. . ,d— 1, then there is a stratified curvilinear 

polyhedral domain uj p C S^ 1 ^ 1 , oJp~ ^ S" 1- ' -1 , and a diffeomorphism (f> p : V p — > 
B d-i x B i guch that = q and 

(17) <f> p {n n Vp) = {rx', < r < 1, x' e uj p } x B l , 

inducing a homeomorphism ft n V p — > {rx 1 , < r < 1, x' € ulj,} X B l of stratified 
spaces that is a diffeomorphism on each stratum. 

The set of singular points of ft then consists of i7(™~ 2 ) and is given as part of 
the definition of ft, but it must contain all the geometric, intrinsic singular points 
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of 9fi. Althought we shall not need this definition here, let us mention nevertheless 
that the desingularization of fi, denoted E(f2), is obtained by gluing in a natural 



way all the sets [0, 1) x uJp~ x B l as in Equation (17). The resulting set E(f2) is 
then a manifold with corners that has a natural structure of a Lie manifold with 
boundary, in the sense of [3]. Then E(f2) — > f2 is a differentiable map that is a 
diffeomorphism outside the set of singular points, in E(f2) the set of singular points 
being the set of points belonging to a face of codimension at least two. 

Let fo(x) > be the distance from x to the set fl^ if x. In general, the function 
fo will not be smooth, we therefore replace tq with an equivalent function r$ that 
is smooth outside QS Q K Therefore, we also have that tq(x) > 0, for x ^ fl^°\ and 
that fo/Vo and tq/tq are bounded functions. We shall say that r$ is the smoothed 
distance to Si' ). We replace then the metric g =: go with g\ :— r^ 2 g. We repeat 
this construction for the remaining non-empty strata in the increasing order of the 
dimension of the strata, each time measuring distances in the new metric. Thus 
rfc is the smoothed distance to Q,^ in the metric gk, and we let gt+i '■= ?~fc 2 <7fc, 
k < d — 2. One can prove that gd-i is a compatible metric on the desingularization 
£(f2) [HH3] and hence we can use the results on Sobolev spaces from those papers. 
Let p := r n ri . . . r^_2- Let us denote by T(£l,TM) the space of restrictions to 
of smooth vector fields on M. The resulting structural Lie algebra of vector fields 
on S(O) is simply V = C 00 (S (fi) ) pT (O, TM) . Thus a basis of V over C°°(S(0)) is 
given by {pdi}. The resulting Sobolev spaces are therefore 

/C?(fi) :={u :P W- a d a u€L 2 (n), \a\<m}^p a - n / 2 H m (n,g d ^), 

where the space H m (Q,h) is the Sobolev space associated to the metric h. Let 
rci(x) denote the distance from x to Q( d ~ 2 K One can prove by induction that rn/ P 
and p/tq are both bounded, so in the above definition of Sobolev spaces we can 
replace p with rn- See [13] for details. 

The fact that the Sobolev spaces £™ are associated to a Lie manifold guarantees 
that Laplacian A satisfies elliptic regularity in the scale of spaces /C™(f2). To this 
end, one also needs to establish that p 2 A — A 9d l is a lower order differential 
operator generated by V and C°°(E(f2)). We also obtain as a byproduct the fact 
that the traces at the boundary of the spaces /C™(0) can also be described in terms 
of the Sobolev spaces on dft associate to the conformally equivalent metric h. 

The Hardy-Poincare inequality can also be interpreted in the setting of the de- 
sigularized metric. Indeed, we have that there exists C > such that every point of 
a; is at a distance < C to the Dirichlet part of the boundary of E(f2) if, and only if, 
there exist no two adjacent faces with Neumann boundary conditions. Then, once 
we know that every point is at a distance < C to the Dirichlet boundary, we can 
prove the Hardy-Poincare inequality in the usual way. 
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